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Abstract. Stochatic excitation of stellar oscillations by turbulent convection is investigated and an expression for 
the power injected into the oscillations by the turbulent convection of the outer layers is derived which takes 
into account excitation through turbulent Reynolds stresses and turbulent entropy fluctuations. This formula- 
tion generalizes results from previous works and is built so as to enable investigations of various possible spatial 
and temporal spectra of stellar turbulent convection. For the Reynolds stress contribution and assuming the 
Kolmogorov spectrum we obtain a similar formulation to those derived by previous authors. The entropy contri- 
bution to excitation is found to originate from the advection of the Eulerian entropy fluctuations by the turbulent 
velocity field. Numerical computations in the solar case in a companion paper indicate that the entropy source 
term is dominant over the Reynold stress contribution to mode excitation, except at high frequencies. 
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1. Introduction 



decade later, entropy turbulent fluctuation has been pro- 
posed as a possible additional exc itation source flStcin fc 
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tion on exeiLaLiun and damping piucesses ul slellai us 

dilation modes. In the solar case, the observed oscillation 
modes are believed to be damped as a result of a com- 
petition between several non-adiabatic and turbulent pro- 



cesses (Osaki, 1990; fioudek et al., 1999). Excitation of so- 
lar oscillation modes is attributed to turbulent convection 
at the surface of the Sun. The first theoretical investiga- 
tion of solar oscillation excitation by turbulent convection 



was by Goldrcich fc Keeley (1977, hereafter GK). These 



authors identified the turbulent term of Reynolds stress 
as the main source term of stochastic excitation of so- 
lar acoustic modes in the wave equation. GK derived an 
approximate estimation for the acoustic power injected 
into the oscillations by turbulent convection which arises 
from an equipartition of energy between the turbulent el- 
ements and the oscillations. The result, however, under- 
estimated the powe r by a factor ~ 10 3 compared to the 



The stochastic mechanism may be understood as fol- 
lows : turbulent motions of the stellar material and tur- 
bulent fluctuations of thermodynamicals quantities, which 
occur in the convection zon e, generate acoustic waves 
( Lighthill, 1952 ; Stein, 1967 ). The corresponding acous- 
tic power then excites resonant modes of the stellar cav- 
ity (oscillations) . This excitation of the oscillations results 
from a forcing by incoherent (turbulent) source terms due 
to the turbulent Reynolds stress and turbulent entropy 
fluctuations. 

An alternative formulation for the power a nd ampli- 
tude oscillation is pr oposed by Balmforth (1992 , hereafter 
B92) and is used by |Houdck ct al. (1999|) to compute am- 
plitude of oscillation for various solar- type stars. There 
is some disagreement between the conclusions of both for- 
mulations concerning the entropy contribution : it is found 



. . . ,i „ \ . — 7777771, . ,.. . . , . . dominant in theoretical investigation (Goldreich et al 
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for sol ar-like oscillating stars have subsequently been com- 

puted by Christcnsen-Dalsgaard fc Frandscn (1983 ) based 
on this simple picture of equipartition of energy between 
turbulence and oscillation. GK considered the adiabatic 
assumption for both the oscillations and the turbulence. A 



1994 , hereaf ter GMK) and as a resu lt of numerical 3D 
simulations (Stein fc Nordlund, 1991), but negligible in 
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some other cases ( Balmforth, 1992j )T 

Both formulations are built following the method de- 
veloped by GK and are based on a simplified descrip- 
tion of the turbulent medium. The spectrum of turbulent 
energy in stellar conditions is, however, not well known. 
The best-known spectrum is, of course, the solar spectrum 



2 



(Espagnet ct al., 1993: Ncsis ct al., 1993| ) and observations 
tell us that a Kolmogorov spectrum does not represent the 
complete turbulent solar spectrum (Ncsis ct al., 1993). 

The purpose of the present work therefore is to inves- 
tigate the effect of several possible models of turbulence 
on the excitation of solar-like oscillation modes and to 
establish their signature in power spectra. To do so, a 
formulation must first be established which includes both 
(kinetic and entropic) contributions in a consistent and 
general way. The main goal is to allow the use of any form 
of the turbulent kinetic energy spectrum, of turbulent en- 
tropy fluctuations and the eddy temporal spectrum. Two 
free parameters are introduced to take into account un- 
certainties in the definitions of the coherence time of the 
turbulent eddies and in the magnitude of the wavenumber 
at which the inertial regime starts. 

Once the general formulation is established, we study 
the specific case of a Gaussian time and Kolmogorov en- 
ergy spectra, which enables us to compare our findings 
with results of previous works. In a companion paper, the 
present formulation is applied to the solar case and sev- 
eral turbulent models found in the literature are discussed. 
Both free parameters can be calibrated with solar obser- 
vations and used to compute oscillation amplitudes for 



othe r potent ial solar- like oscillating stars ( Samadi et al 
2000 |, BOOlbl ) in view of forthcoming space seismological 



experim ents : COROT (|Baglin fe The Corot Team, 1998Q 



MONS QKjeldsen fc Bedding, 1998|), MO ST ( [Matthews 
1998|)7lEDDINGTON ( |Favata et al.,2000| ). 

Sect. 1 recalls how the stochastic excitation mecha- 
nism of stellar oscillations is modeled by source terms in 
the wave equation. We assume the modes are damped and 
find that they are excited through t he turbulent Reynolds 
stresses ( Goldrcich fc Keeley, 1977 ) and the turbulent en- 



tropy fluctuations arising from the turbulent nature of 



the stellar convection region (Stein fc Nordlund, 1991 ; 
Balmforth, 1992| ; jGoldreich et al., 1994| ). We adopt the GK 



procedure and assume that the oscillation modes are de- 
coupled from the turbulent medium. This description gives 
rise to an inhomogeneous wave equation for the oscillation 
modes. In this framework, the acoustic turbulent source 
acts as a forcing term for the oscillation normal modes. 



We consider adiabatic radial oscillations in the Cowling 
(194l[^approximation. A homogeneous, isotropic, station- 
ary turbulence is assumed. A formulation for the Reynolds 
stress contribution which can include any type of turbu- 
lent spectrum is established in Sect. 3. In Sect. 4, the con- 
tribution of turbulent fluctuations of entropy is worked 
out. We find that it is the advection of the entropy turbu- 
lent fluctuation by the turbulent velocity field which yields 
the entropy source. An appendix shows that cross terms 
between Reynolds and entropic sources do not contribute. 

Sect. 5 finally establishes a general formulation which 
can be used to investigate any type of turbulent spec- 
trum. Finally, Sect. 6 discusses the complete formulation, 
its limitations and advantages. 



2. Turbulent stochastic excitation 

2.1. The inhomogeneous wave equation 

In the basic hydrodynamical equations, we use the sym- 
bols P, p, v and g to denote respectively pressure, density, 
velocity and gravitational acceleration. Equilibrium quan- 
tities are represented with a subscript 0. Each variable /, 
except for the velocity v, is written as the sum of the 
equilibrium quantity , fo and a Eulerian fluctuation, f±, 
f = fo + fi and we retain terms which are linear and 
quadratic in the variables Pi and p\ and neglect the grav- 
itational perturbation. Accordingly, one obtains for the 
perturbed momentum and continuity equations: 

dpv 
~dt 



+ V : (pvv) + VPi - Pl g = 



^ + v.( H = o 



(i) 

(2) 



The perturbed equation of state in a Eulerian descrip- 
tion is given by: 



P 



C 2 s Pl 



ot s s\+ a pp p\ J r a ss s\+ oipspisi (3) 

where s is the entropy, a s = (dP/ds) p , c s = T% Po/po de- 
notes the average sound speed, Ti = (91nP/91np) s is the 
adiabatic exponent and a pp , a ss and a ps are the second 
partial derivatives of P versus s and p. 

We assume adiabatic oscillations : the Lagrangian en- 
tropy fluctuations are therefore only due to turbulence. 
Lagrangian and Eulerian entropy fluctuations are related 
to each other by 

dSs t dsi 
~dT ~ ~dt 

where Sst is the turbulent Lagrangian entropy fluctuation. 
The subscript t refers to turbulent quantities. For later 
use and with the help of Eq. (^) , the time derivative of the 
Eulerian entropy fluctuation si is written as 



v.V(s Q + si) 



(4) 



dSst 
: ~dT 



a s v.Vs - V.{a s siv) 
+s\V.'Va s + a s si V.v . 



(5) 



-L)v = 



(6) 



The velocity field v is split into a component due to the 
pulsational displacement Sr osc and a turbulent component 
«asi) = v osc + u. 

Linearisation of Eq.([l]-||) yields for the velocity field, 
in the absence of turbulence (u — 0), the homogeneous 
wave equation 

^dt 2 

with the linear wave operator : 

L(X) = — [V (c 2 s V.(p X) + a s X.Vs ) 
Po 

-g V.( Po X)} . (7) 

With appropriate boundary conditions ( Unno ct al., 1989| ) 
one recovers the usual eigenvalue problem : 

L(€(r)) = - ojI £(r) (8) 
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where luq is the oscillation frequency and £(r) is the adi- 
abatic (real) displacement eigenvector. 

In the presence of turbulence, the pulsational displace- 
ment and velocity are written in terms of the above adi- 
abatic solution £(r,t) and an instantaneous amplitude 
A(t). Accordingly 



^ = i(-i W „#){(r)e-- ( + a; ) 



(9) 



where cc means complex conjugate. 

Differentiating Eq. (Q) with respect to t, neglecting 
non- linear terms in v osc , assuming an incompressible tur- 
bulence (V.it = 0) and using Eq. (||, [| ||) yields the in- 
homogeneous wave equation 



Po 



dt 2 



Po' 



T>{v osc ) = 



with 



ft = -V : {p uu) 

the turbulent Reynolds stress, 



(10) 



(11) 



d d 

■7-Vh t = -V ( a s —8s t - V.(a s usi) + s lU .Va s ) , (12) 
at at 

the source term due to the turbulent entropy fluctua- 
tions. The first term in the RHS of Eq. (12) is due to 
the Lagrangian entropy fluctuation which was considered 
by GMK and B92. The last two terms are due to the buoy- 
ancy force associated with the Eulerian entropy fluctua- 
tions. These terms contribute to the excitation as much as 
the Reynolds source term. The time derivative operator in 
the LHS of Eq. jl^) in the definition of h t is introduced 
for convenience. 

The operator T> involves both the turbulent velocity 
field (it) and the pulsational velocity and is defined as 



d 

V(v osc ) = 2— (V :puv osc ) + V(a 



Pi 



• Vs 



dt \ po 

-V.(a s v osc si) + siV osc .Va s (13) 

+a s si V.v osc - (a ps si + 2a pp pi) V.(pu osc )J . 

As will be shown in the next section, this term contributes 
to the dynamical linear damping. 

We assume a "free turbulence", i.e. the turbulent 
medium evolves freely and is not perturbed by the oscilla- 
tions. The continuity equation is verified by the turbulent 
medium such that 



dp t 
dt 



(14) 



As a consequence of neglecting the oscillation perturba- 
tions in higher order terms, the Eulerian fluctuations Si 
and pi can be seen as due only to the turbulence in evalu- 
ating the quantities V/i t and T>. We therefore substitute 
pi (resp. si) by p t (resp. s t ) in Eq. ||P). 



It can easily be shown that the additional terms ap- 
pearing in the RHS of Eq.([lO|) are of order M 2 , where M t is 
the turbulent Mach number, compared with the Reynolds 
source term (see also GK). As M t is small in the solar 
convection zone (M t < 0.3), these additional terms have 
been neglected. 

The wave operator {d 2 /dt 2 — L) acting on the turbu- 
lent velocity field u in Eq. (|l0|) gives rise to contributions 
which are either negligible compared with the Reynolds 
source term or of the same order. In this last case, the 
associated source term does not contribute to the wave 
excitation because it is linear in terms of the turbulent 
fluctuations (see also Sect. 4.1). Therefore, in the LHS of 
Eq.([l0]), v can be replaced by v osc . 

In deriving Eq.(||), the time variation of the amplitude 
A is neglected since we assume, as in the solar case, that 
the eddy time correlation is of the order of a few minutes in 
the excitation region and the associated angular frequency 
is comparable to too (~ 3 mHz), which is much larger than 
the oscillation damping rate rj (~ 100 pRz). 

2.2. Mean square amplitude 

Substituting Eq. @ into Eq. (|l(i|), gives, with the help of 
Eq. (|), 



2tol^-p Z + iLU T>(A£) 



dt 



+iu> Po 



ith 



dt 2 



El 

Po 



A£ 



as 

dt 



S = -(f t + Vh t 



(15) 



(16) 



As previously mentioned, periods of oscillations are 
much shorter than their lifetime (~ 1 /rj) such that loo 3> rj 
and thus \d\o.A/dt\ ^loo- Therefore the second derivative 
of A arising in Eq. (|l5|) has been discarded. Multiplying 
Eq. ( [l5|) by po£*(r,t) and integrating over the stellar vol- 
ume gives 



dA 
~dt 



Act A = 



1 



2cj2/ 



OS 

' dt 



(17) 



where we have defined the complex quantity Act — iAoj + 
rjD which comes from the contribution of the operator 
2?(£). This expression is linear in the eigenfunction £, 
hence Aui corresponds to a 'dynamical' shift of the oscil- 
lation frequency due to turbulence. This shift is expected 
to be smaller or of the order of non-adiabatic frequency 
shifts, which are not considered here. The dynamical shift, 
Aco, will therefore be neglected. The real part, r/D, con- 
tributes to the dynamical damping. As in B92 and GK, 
rj£, is replaced in Eq. (|lj) by the global damping rate rj in 
order to take into account all damping processes. The so- 
lution of Eq. (0) after integration by parts over the time 
is given by : 

A(t) = % -?— [ dt' f d 3 xe^ +l ^ t ' ^(x).S(x,t') (18) 

2UJ I J-oo J 
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and for the the square amplitude: 



\A\\t) = 



dt x dt 2 \ a I rid I r 2 e' 



<r(ti,t 3 ) 



4(^ /) 2 

x ^(n) . 5(n, ta)) (^(r 2 ) . S*(r 2 , t 2 )) (19) 

with cr(ti, t 2 ) = rj(ti + t 2 ) + iu)a(ti — £ 2 )' Spatial integra- 
tions are performed over the whole stellar envelope. For 
convenience we define the following coordinates 
r 2 +ri ti + t 2 

2 2 (20) 

r = r 2 - ri r = t 2 ~t 1 

Xq and to are the average time-space position where the 
stochastic excitation is integrated whereas r and r are re- 
lated to the local turbulence, t corresponds to the fast 
time scale associated with the eddy lifetime and to to the 
slow time scale associated with the oscillation growth rate 
(rj). In the following, Vo is the large-scale derivative asso- 
ciated with xo, V r is the small-scale one associated with 
r and the derivative operators Vi and V2 are associated 
with ri and r 2 respectively. The mean square amplitude 
can be rewritten in terms of the new coordinates as 

1 ft r 2(t-t ) 

1 dice 2 *""*) / dr 

> J2{to-t) 



(Ml 2 (*)) = 



d 3 x d 3 



(21) 



T 

The operator (.) denotes the statistical average performed 
on an infinite number of independent realizations. 

In the excitation region the eddy lifetime is much 
smaller than the oscillation lifetime (~ I/77) such that the 
integration over r can be extended to infinity. We assume 
a stationary turbulence, therefore the source term S is in- 
variant over time to- Integration over to in Eq.(|2l]) and 
using the definition of S in Eq.jl^) yields : 

(\ A I") = i^na ( C * + °S) ( 22 ) 



with 



r 1 - 

U R — 



d 3 x 



dr 1 



d A r 



(23) 



the turbulent Reynolds stress contribution, 

p-\-oc 



a 



d 3 x 



dr e 



d 6 r 



x<(^V.0 1 (^V.0 2 ) (24) 
the entropy contribution. Subscripts 1 and 2 are the values 
taken at the spatial and temporal positions [xn — § , — §] 
and [xo + §, §] respectively. We adopt the Einstein con- 
vention of summation upon repeated indices. Hereafter all 
time integrations over r are understood to be performed in 
the range ] — 00, +00 [. Note that the crossing term between 
the Reynolds source term and the entropy source term 
does not contribute to the excitation (see Appendix |aJ) . 



3. Reynolds stress contribution 

In the derivation of Eq. ( |23| ) and Eq. (|24|) integrations by 
parts have been performed in order for the gradient to 
act on the eigenfunction instead of turbulent quantities. 
We next suppose that the terms (poV^ J ) and (poV 2 £ J ) 
in Eq.(|23|) do not change on the length scale of the ed- 
dies. This implies that — 0. Validity of this assump- 
tion will be justified a posteriori below. Consequently, the 
Reynolds stress contribution can be expressed as : 



CI 



drd 3 re u-u ; u, u 



(25) 



where the following notations have been used : 



u' = u(x - -M - o) u " = u ( x o + 2> to + 2^ ' ^ 



The Quasi- Normal Approximation ( Lesieur, 1997 , Chap 
VII-2, QNA hereafter) reduces the fourth-order velocity 
correlations as follows : 



Ku'ju'i'u'^ 



(27) 



This approximation remains strictly valid for normally dis- 
tributed fluctuating quantities. As shown by Kraichnan 
|(1957 ) and Stein (1967 ) neglected terms can be large and 
therefore not negligible. Here the QNA approximation is 
nevertheless assumed valid as it is found justified when 
considering 3D simulations of the solar convection zone 
(Samadi, 2000, work in progress). 

We denote (j>ij(k,u>) as the well known Fourier trans- 
form of the second-order velocity correlations {u^u") (e.g. 
[Stein, 1967| ). For a stationary, incompressible, homoge- 
neous and isotropic turbulence, the Fourier transform of 
the velocity correlation has the form ( Batchclor, 1970| ) : 



E{k,w) 
ink 2 



(28) 



where E{k,uj) is the turbulent kinetic energy spectrum 
and Sij is the Kroenecker tensor. Following Stein (1967), 
the velocity energy spectrum E(k, uj) is written as 



E(k,u) =E{k)xk{u) 



(29) 



[Stein (1967 ) and Musielak et al. (1994 ) suggest several 
forms for the frequency factor. The Gaussian function is 
the simplest choice and is defined as 



(30) 



The time spectrum is therefore normalized such that 
([Tcnnckcs k Lumley, 1982[ Chap 8.1) 



(31) 
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One obtains from Eq.(12 



d,T(Pre- Mor (uWf 1 -'"" 



K«>l'<) = (2^) 4 y dcod 3 k 
x [<pu(k,w )<i>j m (k,uj) + (f>i m (k,uo)<f>ji(k,u)] . (32) 



The first term in the RHS of Eq. ( |27j ) when inserted into 
Eq. (|||) gives no contribution. With the help of Eq.(|3l|) 
and Eq.(32), the Reynolds stress contribution can be writ- 
ten as 



C 2 R = it 2 J d 3 x (plVi&Vir) J d 3 k 

E 2 {k) 

X (Tijlm + T lJm i) —j-g — Xfc(Wo) 



(33) 



where 



_ , - kik(\ ( kjk m 

J-ijlm - [Oil p- I I Ojm 



Excitation by convection takes place at the top of the con- 
vection zone. In this region, eigenvectors of acoustic modes 
with high radial order can be considered essentially as ra- 
dial and propagating in a plane-parallel manner. Hence 
Eq. ([33j) becomes (see Appendix |C]) 



r 2 



4tt 3 £ 



M 



dm po 



dr 



F 2 (k\ 
d*-j^X*M(34) 



where Q is an anisotrop ic fa ctor similar to Cough's ( 1977 ) 
one, as defined in Eq. (C.6). Because of the term Xk(^o), 
eddies with u>f. > ujq contribute the most to the integra- 
tion over k, i.e. to mode excitation. This is in agreement 
with Goldreich fe Kumar (199C ), who state that acoustic 
emission arises from eddies with lo^t^ < 1. We have 



ku k 



^0 k osc C s ^osc 



-M t > 1 



(35) 



4. Contribution of entropy fluctuations 

As in Sect. 3 we use V'U 3 - in Eq.(g|) to obtain : 

C 2 S = [ d 3 x (V^) 2 / d T e-**> T I d 3 r{h' t h'l) (36) 



where 

h' t = h t (x - §, - J) and h'l = h t (x + §, +§) . (37) 

Integration over the time of the first term in the RHS of 
Eq. Q gives 



ht{x, t) = -a s Ss t + q t (x, t) 
with 

q t (x,t)= / dt' (V.(a s us t ) - s t u.Va s ) 

J —OO 



(38) 



(39) 



Contribution from the crossing term between the term 
a s Sst (the linear term) and the term q t (the non linear 
term), i.e. (a s Ss' t q' t '), vanishes (cf Appendix ^). We are 
left with 



C% = J d 3 x (V^) 2 J dre-^T J d 3 
x ~((a s Ss' t ) MO) + (q' t q':\ 

4.1. Contribution of the linear term 



(40) 



We consider the contribution of the Lagrangian entropy 
fluctuations (i.e. the first term in the RHS of Eq. (|4C|)). 

In the Boussinesq approximation, the entropy fluctua- 
tion St is related to the temperature fluctuation T t as 



Si 



df) p Tt 



(41) 



In a free turbulent medium the Eulerian temperature fluc- 



tuation acts a s a passive scalar (Tennekes & Lumley, 1982 



Lcsicur, 1997 , Chap VI-10). Thus one can expect that the 



where k osc is the wavenumber of the mode, Uk the veloc- 
ity of an eddy with wavenumber k and M t the turbulent 
Mach number, which is small. Eddies with k 3> k osc are 
those which contribute the most to modal excitation. The 
oscillation and the contributive eddies are then well de- 
coupled. Moreover, the stratification does not affect the 
turbulent emission QGoldreich fc Kumar, 1990| ). The above 
comments altogether justify the assumption- at the begin- 
ning of this section- that (PoVi£ J V§£ ) remains constant 1 
over the length scale of the contributive eddies. These con- 2 
elusions also justify the use of a homogeneous turbulence 
because the stratification and the oscillations have a char- 
acteristic scale length much larger than the contributive 
eddies. 



Eulerian entropy turbulent fluctuation acts as a passive 
scalar. 

Let 4> s (k,uj) be the Fourier transform of the correla- 
tion product of the Eulerian entropy fluctuation. F or any 
passive scalar, one has the relation ( Lesieur, 1997 , Chap 
V-10) 



,(k,uj) 



E a (k,u) 

2nk 2 



(42) 



where the scalar spectrum E s {k, to) is related to the scalar 
variance as (Lesieur, 1997, Chap V-10) 



{s 2 (x ,u)) 



dkE s (k,u) 



(43) 



Because the Lagrangian entropy fluctuation acts as a 
passive scalar as well, Sst scales as st- The contribution 
of the Lagrangian entropy fluctuation is therefore propor- 
tional to 



+00 



dr / d 3 n 



-ibJQT / „U\ 



(44) 



where the following notations have been defined 



s t (x 



. to 



L) 

2 I 



fl*(a;o + 5,tb + 5). (45) 
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As was done with the kinetic energy spectrum, the scalar 
energy spectrum E s (k,w) is decomposed such that 



operators Vi and V2 are replaced by the large-scale gra- 
dient Vo- This finally gives 



E„{k,u) =E s {k) Xk {u) 



(46) 



with the same frequency-dependent factor as in 

Sect. 3. Let Xfe( w ) take the Gaussian form : as u>k oc kuk we 
have Xfc(^o) = for fc = 0. As E s (k ) increases as k 2 in the 
vicinity of k = ( Tennekes, 1975| ; Tennekes fc Lumley 
1982[~^hap 8.6), the quantity s (O,a->o) vanishes and so 
does the linear term contribution. We emphasize that in 
the context of the approximation presented in Eq. (^6|) , no 
assumption has been made regarding the behavior of the 
entropy energy spectrum E s (k). Hence, this result remains 
valid for any scalar spectrum and for the velocity energy 
spectrum. 

This result may also be explained as follows: in term 
of mode excitation, the linear entropy source term acts as 
a second-order correlation product (s's") . Turbulence and 
oscillation are coupled through the phase term e~ tuJ ° T and 
through the turbulent time spectrum which is the 

frequency-dependent component of (s's") in the Fourier 
space. Therefore coupling between turbulence and oscilla- 
tion occurs at frequencies close to the oscillation frequency 
uio and thus involves eddies of wavenumber k > k osc M^ 1 
according to Eq. (|35|). On the other hand, the spatial com- 
ponent of (s's") in the Fourier space favors eddies with the 
largest size (k — ► 0). These two opposite effects clearly are 
incompatible and lead to vanishing of the entropy fluctu- 
ation contribution. 



r 2 

L s 



with 



J d 3 x (a s V .|) 2 9 ij JdTe-^ T J 



d 3 i 



r/2 



dU 



-r/2 



dt i ((s t Ui) 3 (sjUj) 4 ) 



(48) 



V* (ln|a. |)V J (m|a s |) 
-2Vj,(lnK |)V£(ln| V .£ |) 
+V*,(ln| V .£|)V J (ln| V 4 



(49) 



Again, we consider the gradient of the stratification as 
radial and a plane parallel approximation in the excitation 
region. One obtains (see Appendix ^ and Appendix ^|) 



C 2 - 



with 



4tt 3 H 



dk 



,3 ( d & 

X ° \ ~d~ 

E s (k)E(k) 



da; Xk{u + uj)xk{u) (50) 



d In 1 



In 



d ln 1 d_^ 
d\nr dr r 



(51) 



and TL an anisotropic factor similar to Cough's anisotropic 
factor (see Appendix |(] 



This does not happen for Ihe eoiilribulion of Ihe 
Reynolds source Lerni, whieh involves Lhe fuurLh-oi'dei' ve- 

locity correlation product. According to the QNA this 
term can be decomposed in terms of a product of two 
second-order velocity correlations. Coupling with the os- 
cillation then becomes non-linear and leads to an effective 
non-zero contribution. Thus, only non-linear terms can 
contribute to mode excitation while linear terms do not. 
This may be considered as a general result and justifies 
neglect of several source terms in section 2.2. 

4.2. Contribution of the non-linear terms 

As the linear term does not contribute to the acoustic 



The Mixing-Length Theory (hereafter MLT Bohm 
Vitense, 1958^ pox, 1968| ; |Gough, 1977| ) provides an es- 



timate of the vertical velocity of the convective flow. 
The corresponding kinetic energy is transferred to smaller 
scales through the turbulent cascade. The kinetic energy 
spectrum E(k) is normalized as 



1 C°° 1 3 

-(u 2 (xo))^j dkE{k) = -<S>w 2 = -u 2 



(52) 



where w is an estimate for the vertical convective velocity 



emissi qn, Eq.(|40|) becomes, with Eq.( 39 ) 



and $ is a factor introduced by Gough (1977) to take into 
account anisotropy effects (Eq. C.6). uq in the RHS of 
Eq.(|2|) has been introduced for convenience. 

The MLT provides a relation between the temperature 
fluctu ations and the vertical convective velocity w ( |Gough 



1977|) 



C 2 S = J d 3 x (Vo-O 2 / dre-^ r J d 3 r 

,t/2 j—t/2 

x dt 3 dti 

J —oc J —oc 

x { Via s . {(s t u) 3 (s t w) 4 ) .V 2 a s 
- 2Via s . ((s t u) 3 V 2 - (a s s t u) 4 ) 
+ (V 1 .(a s s t u) 3 V 2 .(a s s t u) 4 ) } 



gA5 
2$T 



where 8 = 



2\l/2 



(47) 



d In po 
d\nT 



(53) 



Thus Eq. (|l], II) provide an esti 



mate of the entropy scalar variance 

2 



(*?(*(>)) 



IwJ w 



(54) 



where subscripts 3 and 4 refer to evaluations at positions 
[xo — r/2,t 3 ] and [aso + r/2, 44] respectively. Similar way 
to the Reynolds contribution, we have V\a s ~ y\a s ~ 
VQa s . We use an integration by parts and the derivative 



This enables us to normalize the entropy spectrum E s (k) 



2<*?(*o)) 



f°° 1 
J dkE s {k) = -s 2 



(55) 
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Following GMK we define 1Z as 



-6- 



n = -P^ = — ri 



where a c is the mixing-length parameter defined as usual 
by A = a c H p with A being the mixing-length and H p 
the pressure scale height and I\ is the usual adiabatic 
index. The quantity 1Z 2 roughly measures the ratio of the 
excitation by entropy fluctuations to that by fluctuations 
of Reynolds stresses. 

The entropy contribution can then be written as : 



P1V1 

4ir 3 H / dm p u^TZ 2 T 
Jo 

dkE s {k)E{k) 



with T 2 = A 2 g r . 



dr 



up 
Au 



S2„,2 
l 



+00 



duj Xk(uo + u)Xk(u) (57) 



5. Complete formulation 

With Eq.(||) and Eq.© 
Eq.(p2h is rewritten as 



the oscillation amplitude 



<M| 2 ) = 



2 V (lu I) 2 



dk E(k) 



k 2 



(58) 



where 



_ r E ( k ) 



— Xk{oj Q )+U'R 2 T 2 ^4^ L 



u 
Auj 



-Hex:- 



dw Xki^a + u})xk{u) , 



(59) 



where Q and H are given in Appendix 

Res ults in Samadi et al. (2001a| ) for the solar case 



and in Samadi et al. (2001b| ) for Procyon are based on 



the above general expression for the mean square ampli- 
tude. In order to compare with results of previous works, 
we next consider specific temporal and energy spectra, 
namely the Gaussian time spectrum and the Kolmogorov 
spectrum. 

5.1. Gaussian time spectrum 

Let consider the Gaussian time spectrum given by Eq. ( |3(i| ) . 
This time spectrum corresponds in the time space to a 
gaussian function where linewidth is equal to 2rk and Tk 
is the characteristic time correlation length of an eddy of 
wavenumber k. Hence ujk and Tk are related to each other 



uj k = — ■ 

Tk 

As in B92 we define Tk as 
A 

ku k 



(60) 



(61) 



where the velocity Uk of the eddy with wave number k 



is related to the kinetic energy spectrum E{k) as (Stein 



(56) 1967) 



2/, 



dkE{k) 



(62) 



The factor A in Eq.(61 ) is introduced in order to gauge the 
definitions of Tk and Uk which involve some arbitrariness. 

Let us define the wavenumber ko as the wavenumber 
of the largest eddy in the inertial range. Thus we relate 
kg to the mixing length as follows : 



k 



2tt 



(63) 



where the parameter [3 is introduced here again in order 
to gauge the definition of fco- 

For convenience we define the following variables : 



K = — , u K = — 

fc Uo 



LOA 



2ttuo 
A 



The mean square amplitude of Eq. 
form of Eq. (pCf) is given by 



(64) 

for the Gaussian 



(Ml 2 ) = 



A/3 4 



J dmp u A I — 



32ttt](lu I) 2 

x [gs R {x) + nn 2 T 2 s s (x)\ 



with the "source functions" 



Sr{X) 
Ss{X) 



1 



x- 2 



dK 1 k 2 E 2 {K) _ Z 2 
K 3 u K u% 
dK 1 



2(2tt) 



-</2 



K 3 



UK 



k E s (K) \ ( kpE{K) \ 



(65) 



(66) 



(67) 



where X(m, u>o) = u>o/wa and the ratio Z = \(3X/2Kuk- 
As shown below, the ratios fco Es{K)/s 2 and fco E(K)/uq 
only depend on K. 

The amplitude (Eq. |65| ) is very sensitive to the param- 
eters (3 and A because it scales as A/3 4 and the quantity 
A/3 is involved in the exponential function of Eq.(|66]) and 
Eq.([37|). In section 5_^2 some physical arguments yield a 
crude estimate for the parameter (3. Values of parameters 
A, (3 are discussed in some details in Samadi et al. (2001a ) 
in connection with solar seismic observations. 



5.2. Kolmogorov energy spectrum 

The normalization conditio n Eq.(|52]) allo ws us to express 
the Kolmogorov spectrum (Lesieur, 1997, Chap VI-4) as 



E{K) = ^ K-^ 3 for K > 1 
fc 



(68) 



The expression for the mean square amplitude (Eq. 
) involves also the entropy spectrum E s (k) which must 
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therefore be determined. We recall that the energy spec- 
tra of the temperature and of the entropy fluctuations ex- 
hibit the same behavior (see section 4.1). Turbulence the- 
ory predicts that the energy spectrum of the temperature 
fluctuations follows the same scaling law as that of the 
kinetic energy spectrum E{k) in the inertial-convective 
range (Lesieur, 1997, Chap VI-10) and it decreases as 
fc -17 / 3 in the incrtial-conductive range because molecular- 
conductive effects are predominant. This behavior seems 
to be quite well supported observations of solar granula- 
tion as is claimed by Espagnet et al. (1993) and Nesis et al 
(19930 

We next turn to a specific case where we assume that 
the entropy energy spectrum lies in the inertial-convective 
range i.e. the turbulent entropy spectrum scales as the ki- 
netic spectrum. According to the normalization conditions 
of Eq.(5^) and of Eq.(55), this assumption provides the re- 
lation: 



E s {k,uj) _ 1 E{k,uj) 
P ~ 3 ul ' 

The source functions with Eq. 
( p6| ) and Eq. ([37]) become : 

0.02 



8h and Eq.(B9 



S R {X) 
S S {X) 



K 5 



0.008 
~X r 



dK^e 



K P -(Z/V2) 2 



(69) 
Eq. 

(70) 
(71) 



with X and Z defined below Eq. ( |67| ) . The source function 
Sr(X) (Eq. fn5) is similar to the source function SQ(m, uj) 
obtained by B92 (Eq. 2.20). However the source function 
for the entropy contribution Ss{X) (Eq. fn]) differs from 
the one established by B92 (Eq. 3.9) in that the author ex- 
trapolated GK's formulation of the Reynolds stress contri- 
bution for the contribution due to the Lagrangian entropy 
fluctuations. 

For sake of comparison with the GMK formulation, we 
simplify Eq.(|70|) and Eq.(|7l|) by using the fact that most 
of the stochastic emission occurs from eddies with lj^ > uj 
(i.e. with e~(" / Wfc ) ~ 1). For cjo/^a ^ 1 integration over 
K leads then to the asymptotic forms 



S R oc (u> /wa) 



-15/2 for — > 1 , 

LOA 

UJQ 



(72) 



S s oc (cV^r 19/2 for =± > 1 . (73) 

UJ A 

The asymptotic frequency dependence, Eq. (|72j), is as 
that found by GMK. On the other hand, GMK assume 
the same frequency dependence for the two contributions 
while here the source function Ss (Eq. corresponding 
to the entropy contribution exhibits a steeper slope. This 
result shows us that the power emission from the entropy 
contribution is less efficient at high frequency compared 
to the Reynolds term and thus differs from the results of 
GMK but is consistent with the GK statement that the 
contribution of the entropy fluctuation is larger for the 
long period p-modes than for the short period p-modes. 



5.3. Constraints on free parameters 

We next turn to the parameter /3. The time scale at which 
the convective energy dissipates through the turbulent 
cascade is of order A/w. Thus for stationary turbulence 
an d using Eq.(j52| ), the rate of injection of kinetic energy 
e (Lesieur, 1997, Chap VI-3) can be crudely estimated as 



■Uqw/A . 



(74) 



As E{k) = C K e 2 / 3 fc- 5 / 3 (see [Lesieur, 1997| , Chap VI-4) 
where Ck is the Kolmogorov universal constant, which is 
close to 1.5, and from Eq.(pq) and Eq.(f74[) one obtains 



4 /2 W3)- 1/2 a- 



2^ K 
and then 

/3-y^ 2 ($/3)V 2 



(75) 



(76) 



This suggests that crudely 8 ~ 1.9 with $ = 2, the value 
of the anisotropic factor consistent with BV's MLT. 

The value suggested for Q is somewhat approximate. 
Therefore, as in the case of the parameter A, we consider 
as a free parameter. However, the value of (3\ is con- 
strained by an upper limit. Indeed, let Tk be the cor- 
relation time of the largest eddy in the inertial domain. 
The lifetime of the largest eddies in the inertial range can- 
not be longer than the characteristic time A/w at which 
the convective energy dissipates into the turbulent cas- 
cade. Therefore we must have t^ < A/w and according 
to Eqs. (|Sl| , p2|) evaluated for k = k and Eqs.(|63|,|52|) we 
obtain 0X < 2.7 $ 1 / 2 . 

6. Discussion 

In the present work, a formulation has been established 
which yields the oscillation amplitude of a stellar oscilla- 
tion mode when it is stochastically excited by turbulent 
convection. The main result of this paper is the expression 
for the mean square oscillation amplitude (| A | 2 ) given by 
Eq.([58|) and Eq.([59|). The derivation is based on theoreti- 
cal developments of previous works (GK, GMK, B92) but 
an effort has been made to obtain a sufficiently general 
expression which enables one to implement any type of 
turbulent (kinetic and entropic) spectra and any type of 
temporal spectra for the turbulent eddies. 

For comparison purpose, we next focused on a gaussian 
temporal and Kolmogorov energy spectra; we then ended 
up with the same expression for the Reynolds stress con- 
tribution as obtained in GMK and B92. We must stress 
however that in order to use the same formulation for an 
energy spectrum other than the Kolmogorov one, such as 
for instance the Spiegel spectrum, a general expression 
such as Eq. ([38]) must be used. 

As far as the entropy contribution is concerned, we 
found that the linear term due to the entropy fluctua- 
tion gives no contribution and that it is the advection of 







the Eulerian entropy fluctuations by the turbulent velocity 
field which produces a nonzero acoustic emission. 

In the derivation of the expected mean-square ampli- 
tude (Eq. |58| ), several assumptions and approximations 
have been made. For instance, it has been assumed that 
the oscillations and the stratification are decoupled for the 
eddies which contribute to the stochastic power emission. 
This assumption was shown to be valid and is in agree- 
ment with poldreich fc Kumar (199C ). In addition, we 



have used the plane-parallel approximation, which is valid 
in the excitation region. Other assumptions are based on 
results from Stein's work ( 1967 ), such as the separation of 
the kinetic energy spectrum E(k, uj) in term of a purely 
spatial energy spectrum E{k) and a time- dependent fa ctor 
Xfc(w) for an eddy of wavenumber k. As in Stein (1967 ) the 
QNA has been used. We have used crude approximations 
for estimations of the velocity and the life-time of an eddy 
as proposed by [Stein (1967 ). This led us to introduce the 
free parameter A in the definition of the eddy lifetime. 

The entropy fluctuation has been considered to act as 
a passive scalar and we have extended the separation of 
the kinetic energy spectrum in terms of a purely spatial 
energy spectrum E(k) and a time-dependent factor to the 
entropy energy spectrum. 

The MLT was required in order to estimate the power 
injected in the velocity and entropy turbulent cascade (Eq. 
|65| ). This theory, which assumes the Boussinesq approxi- 
mation, is well known to be a crude approximation. For 
instance the MLT predicts that the characteristic size of 
the largest turbulent element is comparable with the scale 
height of the stratification. This is in contrast with the ho- 
mogeneous hypothesis considered here for the description 
of the turbulent medium. However, it has already been 
stressed that the stochastic emission is not affected by the 
stratification. The use of an homogeneous turbulence is 
therefore valid. 

The size of the largest eddy in the inertial range is es- 
timated by the mixing length (A) according to the MLT. 
We have related the wavenumber of the largest eddy in 
the inertial range (fco) to the mixing length. However, as 
for the eddy lifetime, this relation is rather arbitrary and 
therefore involves uncertainties. We have therefore intro- 
duced an additional free parameter (3. 

We have considered an isotropic turbulence. However, 
effects of anisotropy in the amplitude computaion have 
been partially taken into account. In this way, two 
anisotropic factors H. and Q have been introduced for both 
contributions respectively. These factors have been related 
to 



Gough (1977) anisotropic factor 



It is possible to validate some of these approximations 
by comparing them with results of 3D simulations of the 
solar envelope (Samadi, 2000, work in progress): for in- 
stance, the QNA is found to be reasonably valid. 

To date, and in the solar case, several possible turbu- 
lent spectra can be inves tigated and compared with so- 
lar seismic observations ( Bamadi ct al., 2001a ). But we 
can anticipate that the entropy contribution will be dom- 
inant, as already pointed out by GMK. Comparison with 



solar data allows us to calibrate the free parameters which 
in turn can be used to compute oscillation power spec- 
tra for various solar-like oscillation stars. Indeed, unlike 
the Sun, it is not possible to determine the turbulent 
spectra of other stars from observations of the surface 
granulation. In the prospect of forthcoming space seismic 
projects (COROT, MONS, MOST, EDDINGTON) com- 
parison of theoretical computations with seismic data of 
several solar-like oscillation stars will provide useful con- 
straints on stellar turbulent spectra. 
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Appendix A: Crossing terms 

The entropy and Reynolds source terms involve three 
crossing terms arising from Eq. (|Tj) . These terms are pro- 
portional to 



a r {s u i Ua ) 



(A.I) 
(A.2) 



The QNA gives the following relation 



i i i a u\ 
{u lS u^A 



i i i\ i a u\ 
{UlS ){u i u j ) 



I ll\ I I ll\ 

+{u' l u , ;){s' U '{) . (A.3) 

For an isotropic, homogeneous and incompressible turbu- 
lence we have (s'u'J) = ( [Lesieur, 1997| Chap V-8). Thus 
the three terms in the RHS of Eq. (A.3) vanish. 

The entropy source term introduces a crossing term 
between the linear and the non- linear terms of Eq.(|39j). 
This term is proportional to 



d 3 r (s 1 s"u'l) 



(A.4) 



Third order correlation products of quantities follow- 
ing a normal distribution are expected to be equal to 
zero. Turbulent quantities, such as the entropy fluctua- 
tions (si) and the turbulent velocity field (it), are not 
Gaussian random functions, since such turbulence would 



have no energy transfer between wavenumbers ( Lesieur 



1997). However, for a strongly turbulent medium, as in a 
stellar convection zone, approximation of turbulent quan- 
tities by Gaussian random functions is consistent with the 
QNA and can therefore be considered as valid. 

As a conclusion, the crossing term Crs and the cross- 
ing terms between the linear and non-linear terms of the 
entropy contribution do not contribute to oscillation exci- 
tation. 



Appendix B: Entropy contribution 

The QNA provides the relation 

(?4?4S3S4) = (u l 3 u{)(s 3 S4) + (S 3 UI)(S4UI) 

+(s 3 ui){s 4 u l 3 ) . 



(B.l) 
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The terms of type (UIS3) and (W3S4) vanish (Lesieur, 1997, 
Chap V-8) and Eq. ( |48| ) reduces to 

Cj = J d 3 x (a s V 4) 2 9ij J dre- l ^ T J d 3 r 

dt A (44) <S3S 4 ) • (B.2) 



r/2 



du 



-r/2 



Stati onarity of the turbulence allows us to write 
Eq.([B~2|) as 

/2 

dh 



Cl oc 



dTe- luJ0T d 3 r 



-r/2 



dU (nV(t 4 - *a)) (*t*t(*4 - (B.3) 
which can be expanded in the Fourier space as 
C 2 S oc i 2 *) 3 J dTe- ZbJoT d 3 k 

,r/2 ,-• r/2 r- 

x / dt^ / / duj3 duj4 

J —oc J —oc J 

x fa (fc, W3 ) S (fc, w 4 ) e~ 4 * 3 ^ 3+W4 ^ e "4(^3+^4) _ ( B .4) 



Integrating over £3 and t± in Eq.(B.4) and then over r 
yields : 



Cl oc 



(27T) 4 



d 3 /c / duo <f>ij(k, uj) 4> s ( k, uj a + uj) . (B.5) 



According to Eqs.fl2j, |29|, |42|, ^6|), Eq. ([B.5D can be written 
as 

(27T) 4 



r 2 - 
°s - 



J- y d 3 x (a s V .£) 2 Sij / dw / ^ 



xfSffo-^l EM Xk{LU0+Lo)xk{LO ). (B.6) 



47Tfc 2 V y 



2vrfc 2 



Appendix C: Anisotropy effects 

We establish here the derivation of Eq. (^J) and Eq. ( fsol) 
where the anisotropy factors Q and H. have been intro- 
duced. Thus we evaluate these factors for several cases. Let 
(x, y, z) be the real (anisotropic) coordinates and (x,y, z) 
are isotropic coordinates. We assume an isotropic turbu- 
lence in the horizontal layers while the anisotropy occurs 
in the vertical direction only. The two sets of coordinates 
are then related to each other by the anisotropic factor Q 
as 



x = Qx y — Qy z = z 



(Cl) 



Q is also the ratio of the horizontal to vertical correla- 
tion lengths of the eddies. The equivalent relation in the 
wavenumber space is 



k x — Q k x ky — Q hy k z — k z 



where (k x ,k y ,k z ) and (k x ,k y ,k z ) are anisotropic and 
isotropic wavenumber coordinates respectively. Q — 1 cor- 
responds to an isotropic turbulence. We place ourselves 



in the case Q / 1 but suppose that the energy spectra 
remain unchanged i.e. consistent with an isotropic turbu- 
lence. We therefore model only geometrical effects due to 
non-isotropic shapes of the eddies. 
Evaluation of Eq. p3) gives 



r 2 



4tt j 



dm po 



dj r 
dr 



dk 2 Xk{uo) 
dd sin0 I 1 - 



Q 2 cos 2 



(Q 2 - 1) cos 2 6 + 1 
We define for convenience the factor Q as 
1 / Q 2 cos 2 (? 



Q = j dcos 



1 - 



(Q 2 - 1) cos 2 d + 1 



(C3) 



(C.4) 



This leads to Eq. (§4j). 

As was done with the Reynolds stress contribution , it 
is straightforward to obtain Eq. ( p0| ) from Eq. ( |B.6| ) where 
the anisotropic factor TL is 



H = 



dcos ( 



1 - 



Q 2 cos 2 9 



(C.5) 



(S 2 - 1) cos 2 9 + 1 
|Gough (1977 ) defines the velocity anisotropic factor $ as 



$ = 



(C.6) 



The lifetime of an eddy is the same in any direction ; we 
then have 



QA 



QA 



A 



(Ul) 1 ' 2 (U 2 V Y' 2 



(C7) 



where A is the correlation length of the largest ed dy i n 
the vertical direction. Thus one obtains from Eqs. 
a relation between Q and $ 



(C.8) 



$ = 1 + 2 Q 2 



Integration of Eqs (C.4 , C.5) gives 
1 



Q = 


fl 5/2 




2 


n = 


a 3/2 


with 





2a 3 / 2 - 3^ 



1 



3 atanh(\/a) 



[atanh(v / a) — \fa\ 



a = 1 - Q 2 = (3 - $) /2 



(C.9) 
(CIO) 

(C.ll) 



Table Cl gives values for Q and H for different values of 



(C.2) 
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